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Abstract 

An efficient procedure for generating Gram-Schmidt states on a toroidal surface T 2 is presented. 
As an application of the method, low-lying eigenvalues and wave functions for an electron on T 2 
subjected to a constant magnetic field are determined. 
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Toroidial systems are relevant to fusion heat transfer El 5], v ibrational 0, H, 0, 
electromagnetic |g| and, recently, nanostructure physics |^.TloJll[ H^ . For problems 
restricted to the neighborhood of a toroidal surface, i.e., to those on or near toroidal shells, 
Gram-Schmidt (GS) functions on T 2 could potentially prove useful. As noted elsewhere, 
the calculation of higher order GS states may be time consuming and oftentimes subject to 
numerical error p] when performed in a black-box fashion, particulary when trigonometric 
functions are employed as the primitive basis set and integrations are performed over a non- 
trivial integration measure. In this note we show that there exists a simple algorithm for 
constructing higher order GS functions on T 2 that eliminates these difficulties and proceed 
to apply the functions to the problem of an electron on T 2 in a magnetic field directed along 
the z-axis. 

The geometry adopted here parameterizes a toroidal surface of major radius R and minor 
radius a in terms of cylindrical coordinate unit vectors by 

r(0,0) = W(9)p + a sinflfc (1) 

with W = R + a cos9. Applying d to Eq.(l) gives 

dr = ad6 + Wd^cf) (2) 

with 9 = —sinOp + cos^k. The metric elements can be read off of 

dr ■ dr = a 2 d6 2 + W 2 d(f) 2 (3) 

so that the integration measure is 

^gdq 1 dq 2 -> aWdOdcj). (4) 
It will prove useful in what follows to write the Laplacian derived from Eq. (3); it is 0| 

„, Id 2 sin 9 d 1 d 2 

+ 17T,^3- (5) 



a 2 89 2 aW 89 W 2 , , 

The solutions of the Schrodinger equation derived from Eq. (5) can be taken of the form 
ip n {9)e w ^ with v = integer. A procedure for determining free particle surface wave functions 
ip n {9) has been given elsewhere [3] and extended to several special cases in [3 0, llfij . 
However, if an arbitrary potential V(9, <$) is included in the Hamiltonian, many 9, <fi states 
may be necessary in a basis set expansion to achieve suitable convergence for the eigenvalues 
and wave functions of H(9, 4>). 

The 9^—9 symmetry of the Laplacian allows the solutions of the Schrodinger equation 
to be split into even and odd functions, and the primitive basis set can be taken to possess 
this property, 

u n {9) = —=cos[nO], v n {9) = —=sm[n9]. (6) 

'7T \/7T 



Here for the sake of brevity, we will consider only the even functions in detail and later 
comment briefly on the simple modification needed for the odd functions. Setting a = a/R 
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and F(6) — 1 + a cos#, the simple but key point in what follows is the integrals encountered 
in the GS procedure (irrelevant factors of R and a are dropped) 



are non-zero only for 



t nn = [ u n {d)u n {d)F{6)dd (7) 
Jo 

too = 2, t 01 = a; (8) 

t nn = 1, t n>n+1 = a/2, n > 0. (9) 



Write 



n 



*Pn = ^2 C ^mUm, (10) 
m=0 

and use Eqs. (7) - (9) in the standard GS procedure; the first few GS states are sufficient 
to illustrate the general pattern that emerges {j3 = a/2), 

ipo = uo/V2 = N u (11) 



^1 = N 1 [u 1 -N pN u Q \ (12) 

ij 2 = N 2 [u 2 - N x pN lUl + NtfNiNoPNoUo] (13) 

^3 = A^ 3 [u 3 - N 2 pN 2 u 2 + N 2 f3N 2 Nif3NiUi - N^N^N^N^} (14) 

from which it is apparent that 

cum = (-) n+m N n (N n . 1 (3N n . l )(N n _ 2 (3N n _ 2 )...(N m (3N m ). (15) 

To obtain the c nm the normalization factors must be determined. Consider the k th > 1 
unnormalized state it is easy to establish Nk from 

<$ fc |$ fc >=l-iV fc _i/3JV fc _i (16) 

since all < Ui\u i+ j > vanish for j > 1. Starting then from JVf = (1 — 2 / 3 2 )~ 1 , for > 1 the 
squared normalization for a given $^ is 

= r^vf^ (17) 

In summary, Eq. (17) serves to generate all normalization factors and Eq. (15) the GS 
coefficients. The sinusoidal function coefficients may be obtained by the method employed 
above by starting with c n rather than c 00 and letting Aq = 1. 



3 



It is worth noting that from Eq. (15) the ratio of successive coefficients within a given 
(n, m) series 

C'n.ml C'n.m\- 1 ■N m f3N m (1^) 

allows for truncation approximations to be made for certain problems. To illustrate this 
point, write Eq. (17) in continued fraction form; taking for example A4 (which is the largest 
Nk that can be comfortably typeset), 

Nl = ■ (19) 



1-2/3 2 

Any symbolic algebra program can be utilized to evaluate this expression, but once known 
to a given order of a = 2/5 it need not be evaluated any further. To 0(a 8 ) which occurs at 

A 2 = l + — + — + + + O(a 10 ) + ..., 20 

5 4 8 64 128 V ' K ' 

and a final ^ factor multiplies the expression of Eq. (20) to establish right hand side of 

Eq. (18). Since a < 1, concatenating several Nf3N factors causes the series to converge 

rapidly for larger n and smaller m. For quantum mechanical applications, the kinetic energy 

operator yields an m 2 factor for each term in the expansion given by Eq. (10), making it 

possible to truncate higher states at say, three terms for sufficiently large n (as set by the 

particular problem) 

ip n {9) = c nn cos n6 + c n , n _i cos (n - 1)8 + c n ^ 2 cos (n - 2)9. (21) 

An obvious immediate application for the GS states of relevance to nanoscience 19.11.12] 
is the problem of an electron on T 2 in the constant magnetic field 

B = B k. (22) 



In the Coulomb gauge V ■ A = the vector potential A(#, (f>) = |B x r as expressed in 
surface variables is 

A(M) = ^0. (23) 
The Schrodinger equation (here the spin splitting is neglected) 

1 f h \ 2 

H = ( - V + qA ) \1> = (24) 

2m \ i J 

is more simply expressed if we first define 

7o = B nR 2 

In = — 
1 
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7o 
= — 

In 

2m e Ea 2 
£ = ^~ 
after which Eq. (24) may be put in the form 

( d 2 a sin 9 d a 2 d 2 . 2 d r 2 a 2 F 2 \ 

[^e F 89 + F^d 2 ^ + lT ° a d4> ~ JL ^ + T (25) 

= (H T + e)V = 0. (26) 

The basis expansion functions are taken as per Eqs. (11) - (14) with azimuthal eigen- 
functions appended (the magnetic field term considered here does not cause even (+) and 
odd (-) functions to mix), 




e(M)=^>:4j t: \e-t (27) 

The matrix 

H Tnn = (n\H T \n) (28) 

is then easily constructed since the matrix elements can all be written in closed form, and 
the eigenvalues and eigenvectors determined with a six-state expansion for each ^-parity 0| 
(no truncations were performed). 

Table 1 gives the even function ground and first excited state energy eigenvalues and 
wave functions for several values of v and tq. Table 2 gives the same for the odd functions 
[l8| . Three decimal place accuracy was achieved for the eigenvalues/ wave functions for 
four of the six states generated from H rnn when t = 10, which corresponds to a field of 
2.6 T for a torus with R = 50 nm. Although a large scale treatment of this problem was 
considered outside the scope of this work, it should be noted that because the basis states 
are so simple to generate the only inherent limitation to such a treatment is the matrix 
inversion. Additionally, the matrix elements of the Hamiltonian for an off-axis magnetic 
field that comprises the general case 

B = BJ + B z k (29) 
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can also be done in closed form given sufficient patience 

In conclusion, we have presented an algebraic method to derive GS states on T 2 with 
very little effort. The ability to obtain these functions rapidly may prove of use to problems 
in the areas mentioned in the opening paragraph of this note. As an example of their utility, 
we have used them to calculate the spectra of an electron on the surface of a torus in the 
presence of a magnetic field. 
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TABLE 1: Even function (the (+) superscript has been supressed) ground and first excited state 
energies and wave functions. GS state coefficients not shown in the table are much smaller than 
those given, —v states are trivially obtained from those given below. 



V 


TO 


£o 


# 




*i 











^0 


1.122 


.997^1 - .082^2 + .014^3 





5 


.139 


.840^0 - -533^1 + .099^2 


3.208 


.531^0 + -772^1 - .345^2 





10 


-.747 


.744^0 - .637^1 + .198^2 


7.925 


.466^0 + -289^1 - .751^2 


1 





.249 


-.987^0 - .16201 + .021^2 


1.1663 


-.162-00 + .987^1 + .006V> 2 


1 


5 


1.955 


.882^0 - -469^1 + .052^2 


4.806 


.462^0 + .838^1 - .289^2 


1 


10 


2.429 


.759^0 - -627^1 + .176^2 


10.927 


-.457^0 - -327^1 - .755^2 - -333^ 3 


2 





.795 


.931^0 + -364^1 - .018^2 


3.175 


-.342^0 + .891^1 + .296^2 


2 


5 


4.569 


.973^0 - .221-01 - • 068-02 


7.198 


.217^0 + .973^1 - .063^2 


2 


10 


6.851 


-.803^0 + -587^1 - .105^2 


14.940 


.423^0 + -448^1 - .7747^2 + -247^ 3 



TABLE 2: Odd function (the (-) superscript has been supressed) ground and first excited state 
energies and wave functions. GS state coefficients not shown in the table are much smaller than 
those given. — v states are trivially obtained from those given below. 



V 


TO 


£o 


*o 


£\ 


*i 








.977 


.996^1 - .086^2 


4.033 


-.087^1 - .991^2 + .103^3 





5 


2.510 


.948^1 - .316^2 + .047^3 


5.834 


-.31701 - .914^1 + .250^3 





10 


5.749 


.757V>i - .617^2 + -210^3 


11.128 


.600^1 + .530^2 - -569^3 


1 





1.264 


.999^1 - .037^1 


4.411 


-.036^1 - .997^2 + .067^3 


1 


5 


4.125 


.963^1 - .269^2 + .026^3 


7.441 


-.268^1 - .938^2 + -218^3 


1 


10 


8.760 


.776V>i - .602^2 + .186^3 


14.000 


-.58101 - .567^2 + -5580 3 - .166^4 


2 





2.041 


-.995-01 - .089^2 


5.532 


.088^1 - .995^2 - .046^3 


2 


5 


6.359 


.99101 - .132^2 - .022^3 


9.774 


.12901 + .985^2 - -113-03 


2 


10 


12.710 


-.829-01 + .548^2 - -115^3 


17.611 


.51801 + .676^2 - -513^3 + -107^4 



8 



